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^ ■ Abstract 

l^j ■ We analyze spherical dust collapse with non-vanishing radial pressure, 11, and vanishing tangen- 

C^^ ' tial stresses. Considering a barotropic equation of state, n = 7p, we obtain an analytical solution 

O . 

in closed form — which is exact for 7 — —1,0, and approximate otherwise — near the center of sym- 
metry (where the curvature singularity forms). We study the formation, visibility, and curvature 
VJ ' strength of singularities in the resulting spacetime. We find that visible, Tipler strong singularities 

^r\ can develop from generic initial data. Radial pressure alters the spectrum of possible endstates for 

collapse, increasing the parameter space region that contains no visible singularities, but cannot by 
•1— I . 
j^ ■ itself prevent the formation of visible singularities for sufficiently low values of the energy density. 

H ' 

C^ ' Known results from pressureless dust are recovered in the 7 = limit. 
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1 Introduction 



It has long been known that under a variety of circumstances, spacetimes which are sohitions of 
Einstein's equations with physically reasonable regular initial data, inevitably develop singularities 
|l|. These are events at which Riemannian curvature typically diverges, the spacetime is geodesically 
incomplete, and classical general relativity necessarily breaks down. The question of whether these 
singularities can be visible is one of the outstanding problems in general relativity. 

In an effort to protect the applicability of general relativity, Penrose proposed that such singularities 
might be hidden by an event horizon, and thus invisible to an asymptotic observer, i.e., they cannot be 
globally naked j^]. This constitutes in essence what has become known as the weak cosmic censorship 
conjecture. However, it is quite possible — at least in principle — for an observer to penetrate the event 
horizon and live a rather normal life inside a black hole. This motivated the strong cosmic censorship 
conjecture, which broadly states that timelike singularities cannot occur in nature, i.e., they cannot be 
visible even locally |3|. 

A lack of tools to handle global properties of the Einstein equations (and respective solutions), 
together with their high non-linearity, have been the main obstacle to provable formulations of either 
form of the cosmic censorship conjecture. Whilst efforts are being undertaken in this direction M, one 
can hope that the detailed study of specific models helps to isolate some defining features of singularity 
formation and structure, thereby contributing towards a precise, counter-example-free formulation of 
the conjecture. 

One such model is Lemaitre-Tolman-Bondi (LTB) p inhomogeneous dust collapse, whose general 
solution is analytically obtainable in closed (or parametric) form. The exact solvability of this model 
has resulted in many detailed studies by various authors p, ^1. These analyses show that, from generic 
initial data, a null (central) singularity develops, which is Tipler strong ||], and can be locally or globally 
naked, depending on the differentiability of the initial density profile at the center M . In addition, the 
singularity is stable against initial density perturbations pO[ , and marginally stable against linear non- 



spherical perturbations 1 1 1 . It is fair to say that singularity formation and structure in LTB collapse 
are now very well understood. 

Whilst spherical symmetry may, arguably, constitute a reasonable approximation to realistic collapse 
|l2|, pressureless dust evidently fails to be a realistic form of matter, especially at the late stages of 
stellar collapse, where an effective equation of state must be considered, and radial and tangential 
stresses come into play udk . The inclusion of pressure has been studied analytically in the perfect fluid 



and anisotropic pressure cases; the former, when self-similarity p4| 15 1, or special equations of state 
pql are assumed. 

Early analytical studies of perfect fluid collapse were restricted to approximate solutions of Einstein's 
equations near the singularity p7| ]. This technique was generalized to include a barotropic equation 
of state, p — kp, for perfect fluids, but the results were restricted to the analysis of non-central shells 
due to simplifying assumptions 111]. Numerical studies of self-similar p9| and perfect fluid collapse 
with barotropic equation of state pO| have shown the formation of a central visible (globally naked, 
in the self-similar case) singularity. The work of Joshi and Dwivedi pl| showed the possibility of 
naked singularity formation in the collapse of general "type I" matter fields 123], but their analysis 
does not provide an explicit relation between the initial data — such as densities and pressures — and the 
conditions for formation of visible singularities; (such a relation is very clear in the LTB case). 

Recently, models with only tangential stresses have been analyzed. In particular, the Einstein cluster 
class is now completely understood ]^. Tangential stresses tend to uncover part of the singularity 
spectrum, thereby showing that the role of tangential stresses is not negligible, and lending support to 
the idea that naked singularities develop from generic gravitational collapse ]|2J] . 

In this paper, we present an new solution of spherical dust collapse with non- vanishing radial pres- 
sure, and vanishing tangential stresses. This solution is exact for 7 = — 1 and 7 = 0. For other values 
of 7 G (—1, 1), it asymptotes the exact solution in the r = limit, which is the region of interest, where 
the curvature singularity forms. In particular, the solution can be made to obey Einstein's equations to 
arbitrarily small error, by considering an arbitrarily small neighborhood of the r = 0. In this way, the 
role of radial pressure — as opposed to the combined contribution of radial and tangential stresses, as in 
the perfect fluid case — is isolated and its effects on singularity formation and structure become clear. 
For configurations with radial pressure 11 linearly proportional to the density (11 = 7/)), that can end 
up in either a black hole or naked singularity, an increase in radial pressure leads to a decrease in the 
parameter space area for naked singularities. In the initial two-parameter data space (radial pressure 
vs. density), visible singularities are shown to exist for the entire range of the adiabatic index 7. This 
shows that, whereas radial pressure can cover singularities that would otherwise be visible, it cannot, 
by itself prevent the formation of visible singularities. 

This paper is organized as follows. In Sec. II, Einstein's equations are solved for a spherically 
symmetric system with non-vanishing radial pressure. The ansatz of marginally bound configurations 
yields a complete analytical treatment in closed form. In Sec. Ill, conditions for singularity formation 
are discussed. Section IV studies the visibility of the central curvature singularity and discusses the 



role of radial pressure. Curvature strength is computed in Sec. V. Section VI discusses the pure dust 
limit (7 — 0) of our model. Section VII concludes with a summary and discussion. 
Geometrized units, in which G = c = 1, are used throughout. 

2 Spherically symmetric collapse with radial pressure 

We consider a general spherically symmetric metric, in standard spherical coordinates {t, r, 9, (/>}: 



ds^ = -e2*(*''^)di2 + e-2*(*--)dr2 + R^{t,r)dn^, 



(1) 



where i?(t, r) is the proper area radius and dQ^ = dO^ + sin^ BdcfP' is the canonical metric of the unit 
two- sphere. 

The stress-energy tensor is: 

r„^ = diag(-p,n,E,E) (2) 



where p, 11, and S are the energy density, radial pressure, and tangential stress, respectively. 
With the metric (|l|) the non-vanishing Einstein tensor components are: 
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where ' = dr and ' = dt- 

Introducing the auxiliary functions 
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Einstein's equations can be recast as 



m = AttR^R'p, 
m = -inR^im, 



(10) 
(11) 



$' = ip + ny 
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Note that the function ?7i(t, r) is just the Misner-Sharp mass p5[ : 

m(i,r) = |(l-i?,„i?.b5''^). 



(12) 
(13) 

(14) 



It then follows that the function fc(i, r) is the binding energy per unit mass of a shell r, with area radius 
R(t, r). Gravitationally bound configurations have < A: < 1; unbound configurations have A: < 0, and 
k = corresponds to the marginally bound case. 

From Eq. (0), an "acceleration" equation for the area radius can be easily obtained: 



2* 



TO , , $' 

4^ni? + -^ + (fc-i)- 



R=<PR-e 



When the tangential stresses vanish, T, — 0, Eq. (ha) simplifies to 

$' = -(p + n)-i f2^n + n' 



(15) 



(16) 



2.1 Marginally bound models 

Let us consider the case k — 0, corresponding to marginally bound configurations. From Eq. (g), 
provided R' > (i.e., absence of shell-crossings), it follows that 

R' 



* 



R'' 



which implies, from Eq. (|lj 



$' = 0, 



(17) 



(18) 



i.e., $ = $(i). 

We obtain a closed system by specifying an equation of state, which we take to be of barotropic 
form: 

n = 7p, (19) 

where causality requires 7^ < 1. From Eqs. (ttq) and dla) it follows that 



'^'^'^-R^{t,ry 
where w(i) is an arbitrary real-valued function of r. 
Equations (p^)-(pl]) now read 

to' = AttloR', 
TO = —Attuj^R. 



(20) 



(21) 
(22) 



Equation ( |21| ) readily integrates to 

m{t, r) = AttuR, (23) 

where the "constant" of integration (a function of t alone) was set to zero by demanding regularity at 



the center. Taking the partial time derivative of the above equation and equating it to Eq. (22) gives 

LoR+{l+-i)ujR = Q. (24) 

Which integrates to 

Luit) = luo{t)R-^-\ (25) 

where ujo{i') is an integration "constant" , which can be fixed as uJo{r) — Cr^'^^ , by the choice i?(0, r) — r, 
where C > 0. Hence, 

co{t)^C[-) . (26) 



From Eq. (jlj) we then have 



tC / r \ 1+7 



Note that C characterizes the "strength" of the initial density profile — which is, by construction, positive 
definite, thus enforcing the weak energy condition — whereas 7 G (— l,f) characterizes the relative 
strength of the radial pressure (which can be negative). If we regard 11 as a hydrostatic pressure, then 
y/j is the sound speed, which is complex for negative pressures. Since we do not know how matter 
behaves in high density regimes, such as the strong field regions surrouding the central singularity, we 
allow 7 to be negative, within the limits imposed by the weak energy condition: 

p + n>0^7>-l. (28) 

(Note that causality implies that this inequality does not saturate). 
The mass function becomes 

m(t,r) =47rCi?('-^] ^, (29) 

and the evolution equation simplifies to 

e-2*i?2^87rC(-^)'^\ (30) 

Now, since $ — ^{t), one can trivially rescale the comoving time t to proper time r, via 

dr = e*(*)dt ^ T= /"e*(*)dt + /(r). (31) 

Setting /(r) = 0, Eq. ^) becomes 

R\r,r)=8nC[^y^\ (32) 
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where the dot denotes partial differentiation with respect to r, and R is to be regarded as a function of 
the independent coordinates r and r. Equation (p3) integrates to 

R{T,r) = A^r^ [To{r)-T]^ , (33) 

where 

A = 27rC(3 + 7)2, (34) 

and ro(r) is the proper time for complete collapse of a shell with area radius R(0, r) = r, which is given 

by 

To(r-) = -^. (35) 

Note that Tq > 0, for r > 0, which is a sufficient condition for the absence of shell-crossing singularities. 

2.2 Approximate nature of the solution 

A direct computation shows that Eq. (BS) satisfies the Einstein equations exactly if 7 = 0, or 7 = — 1. 
For other values of 7 e (—1, 1), the Gge component of the field equations does not vanish identically. 
This arises because of the rescaling from t to r, where an unknown function of r was arbitrarily set to 
zero, to render r and r independent coordinates. This lack of exact solvability of the Ggg component 
of the field equations is not detrimental, and can be made precise since 

Gee oc (1 + 7)7r2(i+^)/(3+T)i^(i, r), (36) 

where F{t, 0) oc ^^(5+37)7(3+7) ^{j^qq ^ £ (_x, 1), Ggg is a monotonically increasing function of r, which 
can be made arbitrarily close to zero by considering an arbitrarily small neighborhood of r = 0, which 
is precisely the worldline of the central singularity. Therefore, by restricting ourselves to the vicinity of 
r = 0, Einstein's equations are satisfied to an arbitrarily small error in the region of interest. 

Thus, although not an exact solution for general 7, our solution satisfies Einstein's equations to 
arbitrary precision, and, because it has the correct (exact) limiting behavior for 7 = —1 and 7 = 0, 
it provides a reliable test-bed for the study of the role of radial pressure in gravitational collapse. In 
addition to being an exact solution for 7 — —1,0, Eq. ( p3[ ) captures all the relevant features of well- 
posed dust collapse for 7 G (—1, 1): implosion, satisfaction of the weak energy condition, regular initial 
metric, and trapped-surface-free initial slice. Another example of a test-bed metric that is not an exact 
solution of Einstein's equations, but is useful for its physical content, is the Thorne test-bed metric for 
inspiralling binaries ||26| . Such non-exact solutions are quite valuable in giving analytical, if qualitative, 
insights into physics that otherwise could only be explored numerically. 



3 Conditions for singularity formation and visibility 

The Kretschmann curvature scalar, K. = RabcdR'^'^'^'^ , is 

, , 1 /r\2(i+7) ^ ^ 

/C = 448.^C^^(-) , (37) 

which diverges at i? = 0, thereby signaUng the existence of a curvature singularity at a time tq, where 
R{To,r) ~ 0. Choosing initial data with R{Ti,r) < [note that the choice fc = precludes the t = 
slice from being a moment of time-symmetry, with R{0,r) = 0], a sufficient condition for collapse is 
R{t, r) < 0, for Ti < T < Tq. From Eq. (p2|), we have 

i?(r,r)^--^(l+7)(-) . (38) 

Hence, the condition for complete collapse of a shell r is 

7 > -1, (39) 

which is always satisfied. Hence, provided i?(Ti < 0,r) < 0, implosion always occurs, and the formation 
of a central curvature singularity is inevitable. Physically, our model corresponds to the initial radial 
pressure profile 

n(Ti,r)=7C 



„l+7 



<2g^n(0,r). (40) 

T — Ti 



i^3+7 

As collapse proceeds, H becomes larger (H > 0), but the pressure build-up is never enough to halt 
collapse (R < 0). To this extent, this model is particularly useful for examining central curvature 
singularities — they always form. 

In order to study the collapse of a finite spherical body, we have to introduce a cut-off at some finite 
coordinate radius re, and match the solution along the timelike three-surface Sc : {r = Tc} thus defined 
to a Schwarzschild exterior. Formally, this is achieved by imposing the standard Darmois- Israel junction 
conditions p7| ] to enforce continuity of the metric and extrinsic curvature along Sc. Since our main 
focus is the central singularity and its properties, we shall not go into the details of the matching and 
dynamics of the junction surface; rather, we present qualitative arguments to show that such matching 
can be achieved and that it does not affect the behavior of the solution in the vicinity of the singularity. 
First, we note that the pressure falls off as r^^, and thus we can consider a sufficiently large r^, such 
that n(r, Tc) <C 1 (in appropriate mass-type units). If the cut-off at Tc is sharp (e.g. step- function) , 
one of the components of the extrinsic curvature will fail to be C*^ on Ec- This can be remedied by 
considering a thin shell of matter, r e (re , Tc + e) , where the radial pressure (hence the density) decreases 
smoothly from re to r^ + e, where it vanishes. In this case, the matching with a Schwarzschild exterior 



is trivial at Tc + e, and one only needs to be concerned with the influence of the thin shell on the central 
singularity. This thin shell can be thought of as a perturbation that travels at the local sound speed 
^/7 |l9]. Hence, provided we choose r^ ^ 0{^/jTi), any hydrodynamic perturbations will arrive at the 
center after the singularity has formed. In addition, from a causal viewpoint, the existence of the thin 
shell will only introduce a discontinuity on the AH, when the latter crosses the worldlinc of Vc — the AH 
effectively "jumps out" — collapsing afterwards to the central singularity. 

3.1 Apparent horizon 

The singularity curve is given by Eq. (Bq) . The evolution of the apparent horizon can give insight into 
the causal structure near the singularity. In the adopted spherical coordinates, the apparent horizon 
(AH), which is the outer boundary of a region containing trapped surfaces, is given by 2?7i(Taii(r), r) — 
-R(Tah(r),r). This, from Eqs. (H) and (|3|), gives 

Tah(r) = Mr) [l - e(3+7)/(i+7)] ^ (41) 

with 

e = Vs^, (42) 

where the exponent (3 + 7)/(l + 7) G (2, +00), for —1 < 7 < 1. 

For non-central shells, the mass function m{T, r) is positive definite and finite on any regular surface 
and, from Eq. (p2|), increasing afterwards (m > 0). Therefore, for collapsing configurations, the quantity 
2m/ R increases monotonically, from its initial value (less than unity, which is the condition for non- 
existence of trapped surfaces for regular initial data) , through unity (whence becoming trapped) , before 
diverging positively at the singular surface, r — To{r). The condition for absence of trapped surfaces 
on a given initial slice, 2m/ R < 1 on any Ti < hypersurface, reads 

|^.a.c (!)-<:. m 

Since (r/i?)^+'' < 1, Vtj < and 7 > — 1, the above condition is satisfied provided 

SttC < 1, (44) 

with strict inequality for the 7 = — 1 limiting case. Thus, we have Tshif) < ''0, where the inequality 
saturates for r = 0. Therefore, as expected, non-central shells are strongly censored and the singularity 
cannot be visible even locally. Only the central r = singularity can be naked. Accordingly, we shall 
hereafter discuss the central (r = 0) singularity only. 



4 Visibility 



To show that a singularity is (at least locally) naked, one has to show the existence of non-spacelike 
future-directed outgoing geodesies with their past endpoint at the singularity. In geometric terms, this 
is equivalent to requiring that the area radius increases along such geodesies, i.e., (d-R/d?')oRNG > 0- 
The equation for radial null geodesies (RNGs) is 



RNG 



±R' 
^R 1 



1+7+2 1 



To 



(45) 



r 3 + 7 

where the plus or minus sign corresponds to outgoing or ingoing RNGs, respectively. Along outgoing 

RNGs we have 

^=R' + r(^) =R'il + R). (46) 

dr V^VoRNG 

Using the standard procedure p9|, we introduce the auxiliary variables u, X: 



X 



r" , a > 0, 
R 

u 



In the limit of approach to the singularity we have 



R 



dR 



Xq = lim — = lim 

1 dR 



lim 



(47) 
(48) 



i?^o,r^o ar"^^ dr 

lim ^i?'(l + i?). 

Xq is the value of the tangent to the geodesic on the {u,X} plane, at the singularity. If Xq G 
singularity is at least locally naked, and covered otherwise. Now, from Eq. (BSl) we have 

3 + T- 



T{R,r) 



1 



and 



R'{R,r) 



R 



r-(3 + 7) 
Hence, near the singularity, ro(0), we obtain 

X 



' + ^-(^^^ 



(49) 
, the 

(50) 
(51) 



X, 







lim 



1+7 + 2- 



,(l-a)(3+7)/2 



r^o a(3 + 7) 

^(l-a)(l+7)/2 



X(3+7)/2 



1-e- 



X(l+'^)/2 



(52) 



10 



Clearly, for Xq to be finite at r = we must have a < 1, since 7 e (—1, 1). However, if < a < 1, X 
vanishes identically on the t — surface; hence, we must have a = 1, such that Xq is positive definite 
on the singular surface. 

A self-consistent solution exists if a = 1. In this case, the visibility of the singularity depends on 
the existence of positive roots of the equation 

x„2+7 + i±2ex(^+^)/^ - x(^+^)/^ + e = 0. (53) 

Since — 1 < 7 < 1, this algebraic equation is not a polynomial, and thus solutions cannot in general 
be found analytically. We note, however, the formal similarity with the corresponding equation for 
the general spherically symmetric dust case, where only the lowest-order term appears with a negative 
coefficient ||]. 

Even though the existence of real positive roots cannot be determined for general 7, one can analyt- 
ically examine a sufficiently general class of initial data and draw conclusions from it. Take 7 € [0,1) 
and consider a class of initial data given by arbitrary C [subject to the constraint (p4|)] and 7 = l/{2n), 
where n G N+\{1}. The subset of initial data thus constructed is countable and infinite-dimensional, 
and hence physically significant. Equation ( p3| ) can then be rewritten as 

^4„+i ^ ^!i±le^(6"+i)/2 _ ^(2«+i)/2 + e = 0, (54) 

4n 
Z = Xo^/('"\ (55) 

This is a polynomial of odd-degree and hence it has at least one real root. Now, since 6n + l and 2n+l 
are odd, the second and third terms contain ostensively the square-root of Z. It then follows that the 
real solution(s) must be positive. For this class of initial data there always exists at least one positive 
real root of Eq. (^^, and the central singularity is therefore at least locally naked. An analogous 
analysis for 7 e (—1,0] yields the same results. 

A full numerical evaluation of all the real positive roots of Eq. ( ^3|) is shown in Fig. 1, which covers 
the entire initial data space, {& G (0, 1]} (8) {7 G (—1, 1)}. 

It is instructive to examine the following special cases analytically, as they represent the maximum 
positive, negative, and zero pressure cases, respectively. 

The limiting case 7 = 1 reduces to the cubic polynomial: 

z^ + ez^ - z + e = 0, (56) 

whose solutions are trivial to analyze, for a given Q. Standard methods from polynomial theory show 
that there is a critical value Qmax = 0.300283, above which there are no positive real solutions. For 
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Figure 1: Parameter space for the occurrence of naked singularities. The density increases along the 
vertical axis (O), and radial pressure along the horizontal axis (7). The shaded region corresponds to 
configurations that collapse to a visible singularity. Configurations with initial parameters in the white 
region end up in black holes and contain only spacelike singularities. 



< < 6max, there are three real roots, two of which are positive. This case corresponds to the 
stiffest equation of state allowed. Even in such a case, visible singularities occur. 
The case 7 = — 1 is trivial: 

Xo = 1 - e > 0. (57) 

Hence, the singularity is always visible, since in this case 9 < 1. This is consistent with the notion that 
negative radial pressure can lead to visible singularities that would otherwise be covered [ij] . 

The 7 = case corresponds to a subclass of general inhomogeneous dust models |g^ , and will be 
dealt with separately in Sec. VI. 

There are three relevant regions in the parameter space, defined by critical values of the strength of 
the initial density profile, O (or, equivalently, C): 

Region I: 8 > 0.425343. In this region all the singularities are covered, for positive values of the 
radial pressure. Collapse always ends up in a black hole covering a spacelike singularity, for 7 > and 
& > 0.425343. For negative values of the radial pressure, both black holes and naked singularities can 
form. In the 7 = — 1 limit, collapse always leads to a visible central singularity. The more negative the 
pressure, the larger the number of configurations ending up in naked singularities. 

Region II: 0.425343 < & < 0.300283. For negative pressure, collapse always ends up in a visible 
singularity. For positive pressure, there is a critical value, 0c, below which singularities are visible. 
When the radial pressure is sufficiently large (8 > 8c), the singularities are strongly censored and 
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a black hole is the only final state of collapse. This is in marked contrast with what happens with 
tangential stresses, which tend to uncover the singularity. 

Region III: Q < 0.300283. Collapse always leads to a visible singularity, irrespective of the radial 
pressure. In this case, the magnitude of the radial pressure is bounded from above, and it is not 
enough to censor the singularity. This region shows that radial pressure alone cannot rule out naked 
singularities. 

4.1 Existence of an infinite number of ORNGs 

Each real positive root of Eq. (^) uniquely defines an ORNG with past endpoint at the singularity, 
and a definite tangent value on the {u,X} plane. To examine the possibility of a family of ORNGs 
emanating from the singularity, we briefly describe the method of Joshi and Dwivedi [E9| and refer the 
reader to the original reference for further details. 

Let us consider the ORNG equation on the {u,X} plane: 

'^--(f-x) = '-mX,u)-X], (58) 

du u \ du J u 



where 



Xq = lim - = lim U(X, u) 



Now, let us rewrite the algebraic Eq. ( |53| ) as 

V{Xo) = (Xo - X*)(fo - 1) + f{Xo) = 0, (59) 

where Xq is one of the real positive roots, /o is a constant, and the function /{Xq) contains higher-order 
terms in {Xq — Xq). We have then 

^ = i[{Xo-X*){fQ-l) + S], (60) 

du u 

where S{X,u) = U{X,u) — U{X,0) + h{XQ) vanishes at the singularity. Using the integrating factor 
u^-f«^ Eq. (60) integrates to 



Xq = X* + Buf°-^ + M^o^i / Su-f^du, (61) 

where i? is a constant of integration that labels different geodesies. The last term always vanishes in the 
limit of approach to the singularity (u ^ 0), independently of /q. Clearly, if /o < 1, the second term 
diverges, unless i? = 0, which is required by the existence of the root Xq. In this case, a single ORNG 
departs from the singularity. However, if /o > 1, the second term always vanishes at the singularity, 
irrespective of B, and thus an infinite number of ORNGs (parameterized by B) departs from the 
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Figure 2: Sample plots for roots of the algebraic Eq. (|5^): from top to bottom, for 7 — —0.6,0,0.6, 
respectively. The numerical value of Eq. (|53| ) is plotted on the vertical axes, and Xq along the horizontal 
axes. For a given graph (i.e. 7), each curve corresponds to a choice for 0; one sees that there is a 
critical value Oc, below which there are always two real positive roots. 



singularity. Now, note that if Eq. ( p3[) — and thus V(Xq) — has (at least) two real positive roots, it 
follows that {dV jdXo) = /o — 1 < along one of the roots, and /o — 1 > along the other. Hence, 
in such a case, there is a one-parameter family of ORNG leaving the singularity. A straightforward 
numerical analysis of the roots of Eq. (p3) , shows that for initial data in the shaded region (excluding 
its interface boundary; cf. Fig 1.), there are always two real positive roots (cf. Fig. 2). 

Thus, we conclude that there are infinite number of ORNGs with past endpoint at the singularity. 
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5 Curvature strength 

A defining property of the physical seriousness of a singularity is its curvature strength. A singularity is 
said to be gravitationally strong in the sense of Tipler S if every collapsing volume element is crushed to 
zero at the singularity, and weak otherwise (i.e., if it remains finite). It is generally believed — although 
not yet proven pO| — that spacetime is geodesically incomplete at a strong singularity, but extendible 
through a weak one |^, |3l| . 

A precise characterization of Tipler strong singularities has been given by Clarke and Krolak p^ , 
who proposed (among other conditions) the limiting strong focusing condition: There is at least one 
non-spacelike geodesic, with tangent vector K"^ and afHne parameter A (with A = at the singularity), 
along which the scalar ^P = RabK°-K^ satisfies 

lim A^* > 0. (62) 

A— 

This is a sufficient condition for the singularity to be Tipler strong and corresponds to the vanishing of 
any two- form (for null geodesies), or three- form (for timelike geodesies) defined along such a geodesic, 
at the singularity, due to unbounded curvature growth. 



5.1 Null geodesies 

Let us now consider an ORNG with tangent K"^ = {K'^,K^, 0, 0), where K'^ = ^, and 

K^=^ = R'K"-, (63) 

Jri 

where F can be written as an explicit function of the affine parameter, F = ^(A), obeying the differential 
equation (which follows from the geodesic equation, K°-V aK^ = 0): 

^+F^fA_A_lVo (64) 



In our case, we have 



2to 



R^R 



^^^^ =^ [r) -r^ 



e^X-'-^^, (65) 

r^ 



where Eqs. ([47|)-([48|), and ( |64| ) were used. Hence, 

\2 

lim A^* = e^Xn'^^'' lim 



A^o ° A^o F-2r4 

1 



Xl (66) 
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where THopital's rule was used twice, and [from Eq. (p4)] 

Since Xq > always exists (cf. previous subsection), we conclude that the singularity is Tipler strong 
independently of the details of the initial data. We note, also, that there are infinite number of ORNGs 
with past cndpoint at the singularity, each one of them defined by a unique tangent vector (on the 
{u, R} plane) Xq, which is determined from the initial data (C and 7). This constitutes a rather robust 
result, in that (i) the singularity is Tipler strong irrespective of the initial data, and (ii) it is so along 
an infinite number of ORNGs terminating at the singularity in the past. 

5.2 Timelike geodesies 

Let us consider radial timelike geodesies (RTGs), with tangent vector ^" — ^g^, where A is proper time 
along the geodesic, and 



r = ±VT+R^Wf, (68) 

e-R' + 2CR' + j^iCYR' + ^^3-R" = 0, (69) 

where the first equation is simply ^°^a = — 1, and the second follows from the geodesic equation. By 
inspection, one sees that Eq. ( |69| ) admits the trivial solution 

r = ±1, (70) 

C = 0, (71) 

which leads to 

T = Ti±(A-Ao), (72) 

r = ro = const., (73) 

where the plus or minus sign refers to outgoing or ingoing RTGs, respectively. The outgoing RTG 
departing from the singularity is given by r = 0, and r = tq + A — Aq, and thus does not belong in the 
spacetime. 

The ingoing RTG is given by r = 0, r = tq — A + Aq, where ti = ro(0) = is the time at which the 
RTG arrives at the singularity. For this RTG, we have 
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From Eqs. (p3) and (uT) we get 

R ^ [Ar^+'' {X - Xof]^ , (75) 

which gives 

2 

^ = (76) 

(A-Ao)2(3 + 7)2 ^ ' 

Thus 

m.^(A-A„)^* = ^. (77) 

Hence, the singularity is always Tipler strong along the RTG, irrespective of the initial data. 



6 The dust limit 

When 7 = 0, we have, from Eqs. ( [3^ ) and (p9|): 

i?(r,r)=(^y(ro-r)i, (78) 

where 

m = ■m{r) = AnCr = 47rCp(0, r)r^. (79) 

Equation (^8|) is exactly the same as the one in the marginally bound inhomogeneous dust case ||28|| . 
However, as can be seen from Eq. (|79|), the 7 = limit of our model is a particular case of the 
general inhomogeneous case. The initial mass function cannot be arbitrarily specified, with the only 
free parameter being the constant C. 

At T — To(r), a curvature singularity forms. Events with r > are spacelike and thus covered (cf. 
Sec. HI) . Of potential interest is the r — singularity. The visibility of such a singularity is determined 
by the existence of real positive roots of the algebraic equation (p3). In this case, this equation reduces 
to the quartic: 

P{Y) = Y* + ^Y^ - y + e = 0, (80) 

where Y = ^/Xq. Standard results from polynomial theory can be used to show that there is a maximum 
value Oinax — 0.425343, above which all the roots of the quartic are complex. This result agrees with 
that of Joshi and Singh (JS) for pure dust pq |, upon identification of the mass function in the two 
approaches, whence ones finds [cf. Eqs. (13)-(15) in Ref. p^]: 

©max = 8^G = Ajs = 16/3js = 0.180917. (81) 

For O G (0, 8max], there are two complex conjugate roots and two real positive-definite roots. (When 
= 0, we obtain the trivial solution Y = Xq = and Y ~ Xq — 1, where the former is a degenerate 
case, where the tangent to the ORNG vanishes on the {u,X} plane). 
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As an alternative way of showing the result above, note that Eq. ( pO| ) defines a two-dimensional 
surface that intersects the P ~ plane along a curve 80 (i^), given by solving Eq. (BOh for 8: 

e (Y) 2^^^—^ 

On the P = Q plane, 80 has an absolute maximum at Y^, ~ 0.581029, with 8o(K») = 0.425343. 

From the results of Sec. V, we find that the central singularities in the pure dust limit are Tipler 
strong along an infinite number of ORNGs and at least one IRTG, irrespective of the initial data, in 
agreement with earlier results for general inhomogeneous dust pO| . 

7 Discussion and Conclusions 

The spherical symmetric collapse of dust with non- vanishing radial pressure 11, and vanishing tangential 
stresses was studied. When restricted to marginally bound configurations {k = 0), the field equations 
were analytically integrated in closed form. The fc = ansatz, together with the equation of state 
n — 7p, implies that the initial data consists of two free parameters, 8 and 7, which measure the 
"strength" of the initial density and pressure profiles, respectively. 

As in the pressureless dust case, a central singularity forms, which is Tipler strong along an infinite 
number of ORNGs, and at least one IRTG. Even though, in the timelike case, curvature growth was 
studied along a particular geodesic, using similar causal structure arguments to those of Deshingkar, 
Joshi, and Dwivedi llQ], we can qualitatively show the existence of an infinite number of IRTG with 
future endpoint at the singularity, as follows. The existence of a future-directed timelike geodesic cr(A) 
with future endpoint at the singularity (proved in Sec. V), implies that its chronological past /^[cr] 
is a timelike indecomposable past set (TIP) |33]. In other words, the TIP is generated by a, which 
is timelike and future-inextendible. The locally naked singularity itself — which is the future endpoint 
of cr, cr(Ao) — constitutes a singular TIP which contains the past I^{p) of any point p e a\{a{To)} 
[ p4| . Consider now another point q ^ cr in the past of the singularity and in the chronological past 
of p. Since I^{q) C I^{p), it follows that there exists a timelike curve from p to q, say C, satisfying 
I^[C] C I^[o']. Consider then a small compact ball B — with a suitably defined "radius", e.g., proper 
geodesic distance along cr, d(j(A, Aq) — in the neighborhood of cr(To), partially contained in TIP /^[cr] 
and in the chronological future of p, i.e., Sn /^[cr] = C ^ and I^[C] C I^{p)- For each point x E C, 
there is a timelike geodesic from p to x. Since C is compact, there are infinite number of such points 
that can be joined by timelike geodesies from p G cr. Now, let p be an arbitrary point on a and take 
the limit of approach to the singularity, ^^(A, Aq) ^ 0+; it then follows that there are infinite number 
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of future-inextendible timelike geodesies with future endpoint at the singularity. 

More important than the directional behavior of curvature strength, is the visibility of the singularity. 
We found an infinite number of ORNGs with past endpoint at the singularity, which is therefore at least 
locally naked. The existence of such ORNGs is independent of the initial data. This differs from the 
pressureless innhomogeneous dust case, where the visibility of the central singularity depends crucially 
on differentiability properties of the central energy density distribution, and it is because our model 
reduces to a pressureless dust class with a particular mass function {m oc r). In our model, the initial 
energy density distribution has a single degree of freedom (0), and is C°° everywhere except at the 
center, where it diverges as r^^. However, the mass function is everywhere finite (it vanishes at r = 0), 
and the metric is everywhere regular. We remark that the divergence of the energy density at r = 
does not preclude our model from being a useful test-bed for cosmic censorship. The issue of cosmic 
censorship is the issue of the development of singularities — signalled by the formation of trapped surfaces 
[|l| — from an initial spacelike slice that does not admit trapped surfaces and on which the induced metric 
is regular, in the causal future of that slice. In our model, the initial data (0 and 7) are such that 
the metric is everywhere regular, and there are no trapped surfaces on the initial slice. We have shown 
that, from a large subset of initial data, trapped surfaces — and, consequently, singularities — inevitably 
form in the causal future of such regular initial slices. Accordingly, whereas this model may not provide 
a realistic approximation for the central density, it constitutes a well-defined test-bed spacetime to 
examine the isolated effects of radial pressure on the formation and visibility of the central singularity. 

We find that positive radial pressure tends to supress visible singularities, and negative radial pres- 
sure tends to uncover them. Our negative pressure data space complements the results of Cooperstock 
et al. for non-central singularities |13], to include the r — case. The most significant result is that, 
irrespective of the strength of the radial pressure, there is always a set of non-zero measure in the 
parameter space that leads to visible singularities. Hence, we find that radial pressure alone cannot 
eliminate the naked singularity spectrum. We conjecture that this behavior is not particular to our 
(marginally bound) model, and that it will exist in any spherical system with an equation of state of 
the form H = "/p. In addition, since tangential stresses alone tend to uncover the singularity spectrum, 
and radial pressure by itself fails to fully cover it, we further speculate that when the former and latter 
are taken together into account, visible singularities will still persist. 

In order to confidently establish the role of radial pressure and the combined roles of tangential and 
radial stresses in the final state of gravitational collapse, further analyses need to be undertaken, of 
non-marginally-bound systems, and explicitly including tangential stresses. Efforts in this direction are 
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currently underway p5[ . 
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